SCIENTIFIC 

REPORTS 




SUBJECT AREAS: 

MAGNETIC PROPERTIES 
AND MATERIALS 

SURFACES, INTERFACES AND 
THIN FILMS 

METAMATE RIALS 

TWO-DIMENSIONAL MATERIALS 



Received 
22 January 201 3 



Accepted 
8 March 2013 

Published 
26 March 2013 



Correspondence and 
requests for materials 
should be addressed to 
Z.H.C. (zhcheng@ 
iphy.ac.cn) 



Tuning magnetic anisotropies of Fe films 
on Si ( 1 1 1 ) substrate via direction 
variation of heating current 

Qiong Wu 1 , Wei He 1 , Hao-Liang Liu 1 , Jun Ye 1 - 2 , Xiang-Qun Zhang 1 , Hai-Tao Yang 1 , Zi-Yu Chen 2 
& Zhao-Hua Cheng 1 

1 State Key Laboratory of Magnetism and Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese 
Academy of Sciences, Beijing 100190, China, 2 Department of Physics, Beihang University, Beijing 1 001 91 , China. 

We adopted a novel method to tune the terrace width of Si(l 1 1) substrate by varying the direction of heating 
current. It was observed that the uniaxial magnetic anisotropy (UMA) of Fe films grown on the Si(lll) 
substrate enhanced with decreasing the terrace width and superimposed on the weak six-fold 
magnetocrystalline anisotropy. Furthermore, on the basis of the scanning tunneling microscopy (STM) 
images, self-correlation function calculations confirmed that the UMA was attributed mainly from the 
long-range dipolar interaction between the spins on the surface. Our work opens a new avenue to 
manipulate the magnetic anisotropy of magnetic structures on the stepped substrate by the decoration of its 
atomic steps. 



The structure and magnetic properties of magnetic nanostructures on vicinal substrates have been widely 
investigated for decades 1-5 . The atomically flat terraces separated by steps are generally used as templates for 
preparing various self- organized nanostructures, including regular arrays of nanodots 1 , nanostripes 2,3 , 
atomic wires 4 and ultrathin films 5 . Furthermore, the terraces can be employed to confine magnetic nanodot 
assemblies and form one dimensional (ID) quantum- well states 6 " 10 . For these low- dimensional magnetic materi- 
als, the interplay between quantum confinement and broken symmetry often emerges novel electronic structures 
and magnetic properties. As a result, the magnetic anisotropy of the magnetic nanostructures, which is vital for 
magnetoelectronic technology, can be drastically manipulated by the stepped surfaces 2 " 5 . Recently, we observed 
the tunability in both the magnetic anisotropy and magnetic coupling of Fe nanodots on a curved Cu(lll) 
substrate with varying the terrace width 11 . 

Si(lll) substrate has been often selected to obtain different stepped surfaces 12,13 . Compared to metallic 
substrates, Si substrate is convenient to get a clean surface, and the surface morphology can be manipulated 
by different treating processes 14 " 16 . Moreover, the magnetic nanostructures grown on the atomically flat 
Si(lll) terraces have advantages in spintronic applications as a natural integration with Si-based semi- 
conductor devices. The obstacle resulted from the interdiffusion between the magnetic metals and Si sub- 
strates can be overcome via inserting a buffer layer of copper or iron silicides 1718 . In our previous work 19 " 22 , 
quasi ID magnetic nanodot assemblies were fabricated on vicinal Si(lll) surface with relatively large miscut 
angles (~4°). Unfortunately, the terrace width of these stepped substrates cannot be tuned easily once the 
miscut angle is fixed, which causes difficulty in investigating the effect of terrace width on magnetic prop- 
erties by using identical substrates. Here, we adopted a novel method to tailor the terrace width of Si(lll) 
substrates via controlling the direction of direct heating current passing through the sample. The magnetic 
anisotropy of the corresponding epitaxial Fe films grown on these Si(lll) surfaces was continuously tuned. 
The six-fold magnetocrystalline anisotropy superimposed on a weak UMA can be characterized in the 
sample with the widest terrace, and the UMA will be enhanced with decreasing the terrace width. The 
magnetic anisotropy constants and the origin of UMA were determined via the magnetic hysteresis curves 
simulation and self- correlation function calculation. This work not only advances the understanding of the 
magnetic anisotropy in nanostructures, but also enriches our capacity of fabricating magnetic nanostructures 
and manipulating the magnetic properties for potential applications. 



SCIENTIFIC REPORTS | 3 : 1547 | DOI: 1 0.1 038/srep0 1547 



1 



Results 

Fig. 1(a) illustrates the schematic configuration of the substrate treat- 
ment and preparation process for these Cu/Fe/FeSi 2 /Si(l 1 1) samples. 
During this direct current heating treatment 19 ' 23 , the current passed 
through the substrate, and the direction of current changed from 
parallel to perpendicular to the step. Fig. 1(b) shows a coordinate 
system used in the magneto -optical Kerr effect (MOKE) measure- 
ments and the latter data simulation. The ex-situ magnetization hys- 
teresis curves were determined in the longitudinal MOKE geometry 
at different sample azimuthal angles, (p. The azimuthal angle (p 



(a) 
Cu 25ML 



Fe 20ML^ 




Figure 1 | The schematic configuration of the sample preparation and 
MOKE measurement, and STM images of sample surface, (a) The 

schematic configuration of the substrate treating with controllable heating 
current and film preparation process, (b) The coordinate system used in 
the MOKE measurements and the latter data simulation, (c)-(f) STM 
images of clean surfaces on four Si(lll) substrates I-IV with different 
width of the terrace. The arrows indicate that the direction of the heating 
current changes continuously on each sample, which induce the vibration 
of the width of surface step, (g)-(h) STM images of 20 ML Fe films on 
p( 2 X 2 )/Si (111) of Sample I and IV with the corresponding small scale 
images in the lower right corner. The typical LEED patterns are illustrated 
in the insets of (c) and (g). 



denotes the angle of external magnetic field with respect to the sur- 
face step (Si [1-10]). 

The large scale STM images as indicated in Figs. 1 (c)-(f) dem- 
onstrate the atomically flat Si(l 1 1) surfaces of four samples I-IV with 
different terrace widths. The direction of heating current for each 
sample was denoted by the arrows in these images. The substrate of 
sample I was heated with the current parallel to the surface steps. The 
Si(lll)-7 X 7 reconstructured surface with the terrace width of 
230 ± 30 nm was obtained as shown in Fig. 1 (c). The Si(lll)-7 X 
7 reconstructured surfaces for all the four samples are proved by 
LEED measurements. For samples II-IV, the angles between the 
heating current and the steps vary from 30° to 90°, resulting in a 
decrease of the terrace width. The corresponding terrace widths are 
d n = 100 ± 20 nm, d m = 80 ± 20 nm, d IV = 45 ± 15 nm, respect- 
ively (Figs. 1 (d)-(f)). 

The phenomenon that the variation of terrace width of Si(lll) 
substrate with heating current direction can be explained by the 
surface electromigration of Si adatoms 15 ' 16 . The step bunching due 
to heating by direct current has been well documented on the Si(l 1 1) 
surface with 1° miscut angle 1516 . It was observed that the step bunch- 
ing was influenced by either annealing temperature or the direction 
of current flow. However, the effect of heating current direction 
either parallel or perpendicular to the step on the terrace width is 
not well understood yet. Here, we tuned the step morphologies by 
annealing vicinal Si(Ul) at a fixed temperature while adjusting the 
applied electric field E from parallel to perpendicular to the step. 
When the Si(lll) substrates was treated by direct current heating, 
in addition to the heating temperature, the electric field also influ- 
ences the surface morphology. The Si adatoms would be driven by an 
electromigration force F e =Eq a , where E is the electric field applied 
on the Si adatoms, and g a is the effective charge of Si adatoms. The 
electromigration force could drive the Si adatoms to flow on the 
surface along the direction of electric field. When the direct current 
was applied parallel to the step (0 = 0°), the Si adatoms would flow 
along the edge without changing the step density. Therefore, the 
widest terrace was achieved (Fig. 1(c)). On the other hand, when 
the direct current was applied perpendicular to the step (0 = 90°), 
Si adatoms driven by the largest electromigration force perpendic- 
ular to the terrace edge would diffuse along the electric field and form 
the narrowest terrace width, as shown in Fig. 1 (f). 

In the framework of the surface electromigration model, a scaling 
law was derived 16 



l~N- a (l/F) q 



(1) 



where / is the average terrace width. AT is the number of atoms in one 
step height, and F is the electromigration force perpendicular to the 
steps. Parameters a and q are positive exponents which vary for 
different systems and heating methods. For our samples with single 
atomic steps, Eq. (1) was simplified as 1~{\/F) q . Using the experi- 
mental data of sample II-IV, the value of parameter q = 1.1 was 
obtained. For sample I, when the electromigration force F perpen- 
dicular to the step is smaller than a critical value, the surface terrace 
will be stable and not become wider 16 . 

Figs. 1(g) and (h) display the surface morphology of 20 monolayer 
(ML) iron film deposited on the FeSi 2 p(2 X 2)(1 1 1)/Si(l 1 1) surface 
of samples I and IV, respectively. The epitaxial growth of bcc a-Fe 
film is verified by the LEED pattern (Inset of Fig. 1(g)). Due to the 
small grains of Fe with the surface fluctuation of 2.4 ± 0.3 nm, the 
single atomic steps become hardly visible. From the high resolution 
STM images in the lower right corner of Figs. 1(g) and (h), we can 
observe that the terrace width has no obvious influence on the grain 
size and shape. 

The longitudinal MOKE measurements were performed with 
applied magnetic field at various orientation (p for each sample as 
illustrated by Fig. 1(b). The representative magnetic hysteresis loops 
measured by ex situ MOKE were displayed in Fig. 2. The variation of 
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Figure 2 | Experimental and simulated magnetic hysteresis loops. Typical in-plane magnetic hysteresis loops at different azimuthal angles of 
Sample I-IV ((a)-(d)) characterized by MOKE. Except for Sample I, Sample II-IV ((b)-(d)) exhibited standard easy and hard hysteresis loops at 
cp = 0° and 90°. The corresponding simulated hysteresis loops at hard axis (red lines) are displayed in (b)-(d). 



hysteresis loops with cp indicated in-plane magnetic anisotropy. For 
sample I, the hysteresis loops varied significantly. However no loops 
with the normalized remanence close to 1 or 0 were found at any 
azimuthal angle cp. As shown in Fig. 2 (a), the hysteresis loop with the 
largest and the smallest normalized remanence M r /M s = 0.782, and 
M r /M s = 0.365 appear at cp = 30° and cp = 60°, respectively, indi- 
cating that a six-fold magnetic anisotropy exists in this sample. The 
difference in hysteresis loops for cp = 30° and cp = 90° suggests that 
an additional magnetic anisotropy is superimposed. For the other 
three samples, the rectangular hysteresis loops with M r /M s ~l and 
the sheared loops without remanence were found at cp = 0° and 90°, 
respectively (Figs. 2(b)-(d)). 

For a comprehensive description of the in-plane magnetic aniso- 
tropy of these four samples, the normalized remanence M r /M s as a 
function of azimuthal angle cp of each sample is demonstrated in 
Fig. 3. The complicated symmetry of magnetic anisotropy is observed 
in sample I from Fig. 3 (a). The angular dependence of the normal- 
ized remanence with eight maxima can be well explained by the 
superposition of a six-fold and a UMA contribution. The UMA with 
the easy axis parallel to the steps makes the values of the normalized 
remanence nearby cp = 0° generally larger than that around cp = 90°. 
Moreover, the six-fold magnetocrystalline anisotropy of Fe(l 1 1) film 
modifies the UMA to generate other six peaks at (p = 30°, 90°, 150°, 
210°, 270°, 330°, respectively. As displayed in Figs. 3 (b)-(d), the 
other three samples exhibit similar two -fold symmetric curves of 
the angular dependence of the normalized remanence. The peak 
width of the curves at (p = 0° increases with decreasing the terrace 
width, suggesting that the UMA enhances and consequently conceals 
the weak six-fold magnetocrystalline anisotropy. 

To quantitatively analysis the various magnetic anisotropics, espe- 
cially the UMA, the coherent rotation model is introduced to simu- 
late the hysteresis loops. Besides the UMA energy with the easy axis 
along the steps, the magnetocrystalline anisotropy energy with first- 
and second-order terms, the surface magnetic anisotropy energy, the 
demagnetization energy and the Zeeman energy are taken into 



account. Thus, the total free energy density of the Fe/Si(lll) film 
under an external magnetic field H can be written as 24 

E = — n 0 M s H-fi + Ki (ot\oL 2 2 + ol\ol\ + vl\vl\ ) + K 2 (ol\ ol 2 2 ol 2 3 ) 

( 2 ) 

-K u (ju-h u ) 2 +K d (ju-h d ) 2 -K s (ju-h s ) 2 

where fi is the unit vector of the magnetic moment and M s is the 
saturation magnetization. K\ and K 2 are the first two cubic magne- 
tocrystalline anisotropy constants, with a,- denoting the directional 
cosines of the vector fi with respect to [100], [010], and [001]. K w K s 
and K d represent the anisotropy constants of UMA energy, the sur- 
face magnetic anisotropy energy and the demagnetization energy, 
with the corresponding unit vector of the easy axis in the anisotropy 
h u , h s and n d . It is noteworthy that the directions of h s and h d are 
normal to the Fe(lll) plane and the film plane, respectively. 
Therefore, a small angle exists between h s and h d , which equals 
to the dip angle of surface steps (for the four samples, Pj = 0.07°, f} n 
= 0.17°, p m = 0.21°, P IV = 0.38°). Our modification to Eq. (2) with 
this small angle /? is vital to obtain accurate simulation results in these 
samples with complicated and weak anisotropics. The bulk value M s 
= 1.7 X 10 6 A/m is used, while K x = 4.5 X 10 4 J/m 3 and KJK d = 0.2 
are adopted from our previous ferromagnetic resonance (FMR) ana- 
lysis for the Fe/Si(lll) system 20 . By using K 2 and K u as the fitting 
parameters, various hysteresis loops can be simulated based on 
Eq.(2). The enhancement of UMA with decreasing the terrace width 
can be confirmed by the following fitted results of the other three 
samples: K uII = 5.1 X 10 3 J/m 3 , K uIII = 6.9 X 10 3 J/m 3 , K uIV = 8.8 X 
10 3 J/m 3 . Unfortunately, due to the limit of this simulation based on 
the coherent rotation model, only the simulation results of the 
sheared loops nearby the hard axis can be well fitted with the experi- 
mental loops. Therefore, the accurate values of K 2 and K u in sample I 
cannot be obtained by this method. 

In order to further explore the contribution of the surface mor- 
phology to the UMA, the surface morphology induced in-plane 
magnetic shape anisotropy is estimated by means of the self-correla- 
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Figure 3 | The angular dependence of the normalized remanence M r /M s obtained from the experimental hysteresis loops. Sample I shows a UMA 
superimposed on a weak six-fold anisotropy, and Sample II-IV display a gradually enhanced UMA with the six-fold anisotropy almost concealed. 

tion function method 19 . According to this method, the UMA origi- where S is the film surface adopted in the calculation, and s(r) is the 

nated from the long-range dipolar interactions can be calculated by height distribution of the surface, which can be obtained from the 

using the surface profile data extracted from the STM images. STM images. Fig. 4 (a) shows the spatial distribution of^ (r) calcu- 

The self-correlation function g^f) described the surface rough- lated from the STM image. The data ofgiffl in the plane of the film 

ness as 25 surface display obviously striped structure along the direction of the 




Figure 4 | The typical spatial and in-plane distribution of the self- correlation function gij(r). The spatial distribution of gij(r) with the corresponding 
calculated STM image (a), and in-plane arrangement (b) of g i} ;(r). The value ofgy(r) in the plane of the film surface display obviously striped structure 
along the direction of the step. 
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Table I | The uniaxial magnetic anisotropy K u fitted based on the 
experimental and simulated hysteresis loops, and the magnetic 
shape anisotropy /C s h ape calculated from the surface morphology 
for every sample 



Samples 


1 


II 


III 


IV 


Ku(10 3 J/m 3 ) 
/C shap e (10 3 J/m 3 ) 


5.2 


5.1 
7.1 


6.9 
9.3 


8.8 
12 



where d is the thickness of the magnetic film, and the constant of 
magnetic shape anisotropy induced by the surface morphology can 
be directly given by [AT] 25 , 

K sh ape = ^ fi 0 M 2 s (N xx - N yy ) (5) 

Consequently, the calculated results for the four samples yield the 
following magnetic shape anisotropy constants: K s h apeI — 5.2 X 
10 3 J/m 3 , K shapeII = 7.1 X 10 3 J/m 3 , K shapeIII = 9.3 X 10 3 J/m 3 and 
KshapelV = 1.2 X 10 4 J/m 3 . 

Discussion 

STM images exhibited that the terrace width of stepped Si(lll) 
substrate can be continuously tuned by the direction of heating cur- 
rent. The MOKE magnetic hysteresis loops revealed that the terrace 
width has a significant influence on magnetic anisotropy of the epi- 
taxial Fe films grown on these Si(lll) substrate. The self- correlation 
function calculations demonstrated that the magnetic anisotropy 
resulted from the surface morphology is in good agreement with 
the experimental one (as summarized in Table I). The good agree- 
ment proves that the tunable UMA in this Fe/Si(lll) system is 
mainly originated from the long-range dipolar interaction between 
the spins on the surface with different terrace widths. 

Our work suggests that the magnetic anisotropy of Fe film on 
Si(lll) substrate can be flexibly tuned by direct heating current. 
To our knowledge, the tunability in magnetic anisotropics of Fe 
ultrathin films on Si(l 1 1) substrate via direction variation of heating 
current was not commonly seen in literature. This finding opens a 
new avenue to manipulate the magnetic anisotropy via decorating 
the atomic step of substrates, which can enriches our capacity of 
fabricating magnetic nanostructures and manipulating the magnetic 
properties for potential applications. 

Methods 

The samples were deposited by molecular-beam epitaxy in an ultrahigh vacuum 
(UHV) system with a base pressure of 2. OX 10" 10 mbar. We used the phosphorus- 
doped Si(l 11) substrate with a small miscut angle of 0. 1 0 towards the [11-2] direction. 
The substrates were cut from the identical Si wafer along different orientations to 
control the angle 0 between the atomic step and the heating current in the following 
substrate treating procedure. After the substrates were introduced into the UHV 
chamber, they were degassed at 600°C , followed by a standard direct current flashing 
procedure at 1250°C to obtain the atomically flat Si(lll)-7 X 7 reconstructured 
surfaces with different terrace widths 19 ' 23 . During this flashing process, the heating 
current was directly passed through the substrate. The Fe films were grown using an 
99.999% pure iron wire heated through e-beam evaporation at room temperature 
with a deposition rate of 1.5 ML(monolayer)/min 20 . In order to avoid the Fe-Si 
intermixing, an atomically flat ultrathin p(2 X 2) iron silicide template was first 
prepared through 1.5 monolayer (ML) of iron annealed for 10 min at 400 °C. The 
deposition of 20 ML Fe was subsequently performed, at which thickness the step of 
iron films was still visible. The evolution of surface morphology and structure of the 
samples were characterized by scanning tunneling microscope (STM) and low energy 
electron diffraction (LEED) equipped in this UHV system. Before ex situ 
measurements, 25 ML of Cu was capped to avoid the samples oxidization, and the in 
situ surface magneto -optical Kerr effect (SMOKE) measurements confirmed the 
magnetic properties did not change due to the capping layer. The ex-situ MOKE 
measurements were carried out to investigate the in-plane magnetic anisotropy of the 
samples. 
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